Photonic Josephson effect, phase transitions, and chaos in optomechanical systems 
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A photonic analog of the Josephson effect is analyzed for a system formed by a partly transparent 
mechanical membrane dividing an optical cavity into two halves. Photons tunneling between the two 
sub-cavities constitute the coherent Jospehson current. The force acting upon the membrane due to 
the light pressure induces a nonlinearity which results in a rich dynamical structure. For example, 
contrary to standard bosonic Josephson systems, we encounter chaos. By means of a mean-field 
approach we identify the various regimes and corresponding phase diagram. At the short time scale, 
chaos is demonstrated to prevent regular self-trapping, while for longer times a dissipation induced 
self-trapping effect is possible. 
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Introduction. - The Jospehson effect describes macro- 
scopic tunneling made possible due to inherent coher- 
ence of the underlying quantum many-body state. The 
resulting Josephson oscillations, however first discussed 
for weakly coupled superconductors [![, has turned out 
to be relevant for a wide range of physical systems like, 
Bose-Einstein condensates in double- well traps [2j , spinor 
condensates Q, graphene and Hall systems and in 
nonlinear optics [6j . Lately, the possibility to observe the 
Josephson effect for photons in systems of coupled high Q 
optical cavities has been discussed @, H| . Together with 
related schemes containing arrays of coupled cavities [§] 
or single cavities filled with a Kerr medium [10;], these 
proposals pave the way for genuine quantum many-body 
models to be studied by means of photons. 

Nonlinearity is a basic ingredient for Josephson oscilla- 
tions; in Ref. [7| a three cavity setup is investigated and 
nonlinearity is accomplished letting the photons interact 
dispersively with matter in the middle tunneling cavity, 
while in [8| only two cavities are coupled and the non- 
linearity stems from a gas of ultracold atoms trapped in 
each cavity. In this work we consider the photonic analog 
of Josephson oscillations in an optomechanical cavity sys- 
tem that has been thoroughly analyzed experimentally in 
the recent past 
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12|. Nonlinearity in optomechanical 
systems is an outcome of the interaction between photons 
and the mechanical mirror/membrane, and is thereby dif- 
ferent from the above schemes which utilize atoms or 
Kerr mediums in order to induce the effective interaction 
between photons. Nowadays, the mechanical oscillator 
can be cooled down close to its quantum ground state [l3[ 
making these systems promising alternative candidates 

In our 
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for observing many-body quantum effects 
version of the photonic Josephson effect, the mechanical 
oscillator adds an additional degree of freedom to the dy- 
namics driven by the photon modes, and hence opens up 
for novel phenomena absent in more regular Josephson 



systems. We will indeed demonstrate that the system at 
hand possesses a plethora of phenomena like self-trapping 
(i.e. localization of photons), bistability, and dynamical 
phase transitions between chaotic and regular dynamics. 
We note that the back-action, here induced by the in- 
teraction between the membrane and the photon modes, 
that occurs when the extra degrees of freedom are con- 
sidered has not been taken into account in 0, @] and one 
expects that those analyses would therefore not predict 
structures as chaos and attractors. 




Figure 1: Model system, a partly reflective membrane divides 
a Fabry-Perot cavity into two parts a and b. With a rate 
g, photons can tunnel between the two parts. At the same 
time, the membrane is displaced from its equilibrium position 
(x = 0) due to the radiation pressure imposed by the photons 
scattered upon it, and thereby the two effective photon fre- 
quencies uj a and are simultaneously shifted. Photon losses 
through the two end mirrors are marked by «, and phonon 
decay of the membrane by 7. 



Model system and effective equations. - The cavity 
setup is depicted in Fig. [1] A thin membrane separates 
an optical resonator into two halves, where photons can 
tunnel from one side to the other [ll], [l2[ • The optical 
force acting on the membrane induces a nonlinearity; the 
position of the membrane determines the photon frequen- 
cies which directly affect the Josephson oscillations. By 
considering an optomechanical system with two fixed end 
mirrors, the finesse of the cavity can be greatly enhanced 
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in comparison to a two mirror cavity with one fixed and 
one movable end mirror. In the regime where the pho- 
ton modes follow adiabatically the membrane displace- 
ment 15], non-zero transmitivity of the middle mirror 
brings out constant shifts of the two photon frequencies 
as well as optomechanical couplings of higher order in the 
membrane displacement than the linear one which would 
be the case of a perfectly reflecting membrane [3]. We 
will assume small membrane displacements, i.e. much 
smaller than the photon wavelengths, and the mirror to 
be positioned to render a linear coupling in x [12] • 

In the following we work with dimensionless vari- 
ables, letting the membrane excitation energy, fvjj m , set 
a characteristic energy scale, i.e. giving a length scale 
I = y/ h/niLd m with to the membrane mass. Denoting 
by x and p the dimensionless position and momentum of 
the membrane, a and b (a* and b') annihilation (creation) 
photon operators of the two cavity modes, the Hamilto- 
nian is 1121 



~2 ~ 2 
p X 

Horn = y + y +CJ (h a + fl b ) 



+g (^b + tfa^J + Xx (n a — n&). 



(1) 



Here, we have taken the equilibrium position of the mem- 
brane in the absence of photons (x) = and assumed the 
effective frequencies of the two photon modes to be the 
same and equal to u>, g is the effective tunneling coeffi- 
cient of photons between the left and right regions, A the 
effective optomechanical coupling strength, and n a = a^a 
and hb = wb are the photon number operators of the two 
modes. The number of photons N p h = h a + hj, is con- 
served and we go to an interaction picture rotating with 
respect to LoN p h- In the new basis, the third term of 
H om is lacking and we are left with only two effective 
parameters g and A. 

In the decoupled case, A = 0, and assuming that ini- 
tially all photons resides in one of the modes, as time 
progresses the photons will tunnel through the membrane 
showing perfect Rabi oscillations with a frequency 51 = g. 
This is irrespective of the number of photons, and hence 
after half a Rabi period all photons have been swapped 
to the opposite mode. For A ^ the dynamics become 
much more complex. The shift of the membrane is pro- 
portional to the imbalance of population between the two 
modes, and thereby the dynamics is highly nonlinear. As 
the displacement induces an effective detuning between 
the two modes, for large displacements the tunneling of 
photons will be greatly reduced, i.e. self-trapped. 

We analyze the dynamics within a mean-held approach 
that has proved to accurately describe the experimental 
Josephson effect in Bose-Einstein condensates 0- We 
write down the Heisenberg equations of motion for the 
photon inversion z — (n a — fit,) /iVo, the phase difference 
= arg (o^b) i the photon loss fraction q = (n a + hb)/No, 
and p and x, and treat them as c-numbers to obtain a 



closed set of coupled equations. N is the number of 
photons at time t — 0. This method allows for losses 
of both the photon and phonon helds to be taken into 
account, and we therefore have introduced the quantity 
q. The mean-held equations of motion read 
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x = P - ~K X ' 

1 

P = -x - -p - XN z, 

z = 2g\J q 2 — z 2 sin0 — kz, 

(f) = 2\x — 2g — rna , 

Jq 2 -z 2 

q = - Kq + 2 K ^, 



(2) 



where n is the photon decay rate and 7 the phonon decay 
rate. The photon reservoir temperatures are determined 
by the average number of thermal photons Nth at the fre- 
quency of the two cavity modes, while the phonon reser- 
voir temperature does not enter since we have truncated 
the equations of motions to contain only x and p and not 
any higher orders. 
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Figure 2: (Color online) Time evolution of the mean-field 
photon inversion z(t) (red curves) and membrane position x(t) 
(blue curves). In (a), the system is closed/conservative, 7 = 
k = 0, and the dynamics is chaotic with the inversion making 
sudden jumps as the membrane approaches x — 0. In (b), 
damping of the membrane has been included, 7 = 0.01, « = 0, 
and the system reaches a self-trapped state (an attractor). 
In (c), 7 = 0.01, k = 0.02, and Nth = 0, and as photons 
leak out of the cavity the membrane's displacement vanishes. 
Finally, in (d) the situation is the same as in (c) but with 
Nth = 200 causing the chaotic behavior of detuned Josephson 
oscillations to last throughout the evolution. In all four plots, 
the dimensionless parameters are g = 0.2, A = 0.1, iVo = 1000, 
and the initial conditions xq — po = <f>o — 0, zq = 0.95, 
and q = 1. Thus, the parameters are such that the system 
starts out in the chaotic regime in all examples . Other initial 
conditions show similar behaviors provided the system starts 
out far from any fixed point or attractor. 

Self- Trapping, bifurcation, and chaos. - Several gen- 
eral conclusions can be readily drawn from the set of 
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mean- field equations ([2]); the first two equations describe 
a damped oscillator "driven" by the term proportional 
to z, the third and forth equations represent a pendulum 
with varying length and "driven" by 2Xx. In the limit 
when the shift of the membrane can be approximated 
as proportional to z, one recovers the standard mean- 
field equations of jnotion for a Bose-Einstein condensate 
in a double- well |2j. However, such an approximation is 
in general not justified and the set of equation ([2]) turn 
out to possess a dynamical richness not found in regu- 
lar Bose-Einstein condensate double-wells. The various 
types of dynamical behavior can be divided in: 

(1) Lossless case, k = 7 = 0. In the lossless case the 
system is conservative in the sense that a volume- element 
in phase space preserves its volume in time. The solu- 
tions evolve on four dimensional tori, which become more 
and more deformed the stronger the nonlinearity gets in 
accordance with KAM-theory [l7|. At C = g/X 2 N = 1 
there is a pitchfork bifurcation characterizing a second 
order dynamical phase transition [3]. For C > 1, there 
is one center fixed point, z = x = p = and 4> = (2rt+l)7r 
(n an integer), and the system shows Josephson like dy- 
namics, i.e. far from the fixed point the photons oscillate 
between the two modes. For C < 1, when nonlinear- 
ity becomes important, the above fixed point becomes 
unstable and two additional center fixed points appear, 
z = ±Vl - C 2 , x = -XN a z, p = 0, 4> = (2n + 1)tt. In 
this regime, self-trapping is expected, but is here pre- 
vented by chaotic dynamics, as demonstrated in Fig. [2] 
(a) displaying time evolution of z and x }19j |. The mem- 
brane shows an oscillatory behavior; when it is displaced 
far from x = the effective detuning between the two 
photon modes is large and z shows detuned Josephson 
oscillations, when the membrane approach x = the de- 
tuning vanishes and the photon inversion reveals sudden 
jumps. Such jumps inhibit the system from becoming 
self-trapped, and moreover, small fluctuations in system 
parameters or dynamical variables render completely di- 
vergent long time evolutions as is characteristic for chaos. 
We should point out that there is as well an additional 
fixed point for (j) = 2nir and x = p = z = 0. 

(2) Membrane losses, 7 ^ 0, k = 0. Contrary to the 
above case, for losses of the membrane the system is dis- 
sipative, i.e. a volume element in phase space decreases 
with time [l8| . A similar pitchfork bifurcation as above is 
found at C(l + 7 2 /4) = 1, but this time the fixed points 
are stable. Demonstrated in Fig.[2](b), far from the fixed 
points the system evolution is in general chaotic to start 
with, but after some time it approaches attractors that 
encompass the fixed points for which the membrane is 
locked and z and sin(</>) (we consider sin(</>) instead of <j> 
since the system is periodic in </>) evolve on an ellipse-like 
surface. Physically, the nonzero inversion implies a dis- 
placement of the membrane, while a nonzero loss rate 7 
tend to restore the membrane to its equilibrium position 
x = 0. Since the photon number is conserved, the mem- 



brane attains a state of self-lasing where the light force 
and the friction balance one another. In this respect 
the system becomes self-trapped. We note that similar 
attractors were found for a single mode optomechanical 
system in the linear regime [20|. As long as we do not 
drive the cavity modes, the asymptotic self-trapped state 
is a priori hard to predict because of the chaotic behavior 
preceding the system's convergence of the attractor. 

(3) Cavity and mirror losses at zero temperature, 7 7^ 
0, k ^ 0, Nth — 0. When the system is not externally 
driven and all boson modes possess losses, the only fixed 
point is the one of all modes empty. This stationary 



solution is reached on a timescale T > max{7 



L }- 



Nevertheless, for times t <^T the dynamics will be deter- 
mined by the size of C(l + 7 2 /4) as shown in Fig. [2] (c). 
The decay rate n has been taken as k = 0.1a in agree- 
ment with the experimentally relevant work [12j |. Up till 
t ~ 50 — 100, the dynamics is similar to that displayed 
in (a), i.e. detuned Josephson oscillations interrupted by 
sudden jumps. In this example < 7 implying that 
after a period of chaotic evolution, the larger decay of the 
photon modes causes them to decay into vacuum, while 
the membrane sustains its oscillations for some longer 
times. 

(4) Cavity and mirror losses at non-zero temperature, 
7 / 0, k / 0, Nth 0- When all modes couple to 
their respective reservoirs, an energy source is required 
in order to prevent decay into the trivial vacuum steady 
state solution. This could either be non-zero temperature 
baths or direct pumping Q- For Nth 7^ 0, the equations 
exhibit a nonvanishing fix point for the photon loss frac- 
Also this case includes a pitchfork 



tion q = q Q = 2^ 
bifurcation. When D 



2<; 



> 1 the bifurcation is sub- 



critical; an unstable fix point for C = C(l+7 2 /4)/go < 1 
splits into two unstable and one stable fix point at C = 1. 
When D < 1, the bifurcation is supercritical and one 
stable fix point for C > 1 splits into two stable and one 
unstable fix point at C = 1. For this situation, we find 



as well two saddle node bifurcations at C 



2D 



< 1, 



where the two stable fix points in the pitchfork bifurca- 
tion meet two outer unstable fix points. Using the proper 
physical parameters k and 7 of the examples presented 
in Fig. [2] (c) and (d), the parameter D = 1 is exactly at 
the separation between the two parameter regimes. The 
chaotic behavior we study manifest itself, however, far 
from the critical points and the actual value of D is of 
minor importance in comparison to C ^ 1. When the 
coupling between the membrane and the cavity modes 
are non-zero, the evolution stays chaotic comprising de- 
tuned Josephson oscillations between sudden jumps in 
the inversion, see Fig. [2] (d). We have also studied the 
case of Nth = but with pumping of the two photon 
modes. The dynamics turn out very similar to that en- 
countered for a non-zero temperature bath. 

As is well known, a fixed point analysis does not reveal 
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the full characteristics of dynamical systems 18[ . We find 
that especially for N t h 7^ 0, the long-time solutions of ([2]) 
exhibit a rich structure. As the nonlinearity is increased, 
the long-time solutions cover a series of periodic dou- 
blings in terms of Hopf bifurcations 18[ , to finally enter 
into a regime of chaos. In a driven system, this feature 
should be captured by varying the pump strength, where 
each Hopf bifurcation render a dynamical second order 
phase transition. 

Concluding remarks. - To summarize, we have pro- 
posed the use of optomechanical systems for the study 
of a photonic counterpart of the Josephson effects. The 
system was shown, on a mean-field level, to display 
bifurcations and dynamical phase transitions, as well 
as a chaotic behavior. In this regime, a new phe- 
nomenon of reservoir assisted self-trapping was estab- 
lished. Throughout we used dimensionless parameters. 
Considering exp erimentally relevant parameters for opti- 
cal cavities OE^i one nas typically 10~ 6 < A < 1CT 1 
and 0.2 < g < 100. In this work we choose the upper and 
lower limiting cases of g and A respectively to achieve 
the strongest possible nonlinearity, giving a characteris- 
tic time scale of one fis. However, similar results are 
obtainable with other values on the cost of longer time 
scales. In real units and using common experimental pa- 
rameters [HI, , one finds that the membrane displace- 
ment is several orders of magnitude smaller than typical 
photon wavelengths implying that the linear coupling as- 
sumption is justified. The cavity decay rate k = 0.02 was 
taken from Ref. [13] and should be of experimental rel- 
evance. Utilizing these parameters, Fig. [2] demonstrated 
establishment of both Josephson oscillations as well as 
the chaotic jumps before the role of dissipation becomes 
too important. It is worth noting that for typical cavity 
frequencies, Nth = 200 as in Fig. [2] (d) implies extremely 
hot photon baths, and it is thereby more likely that exter- 
nal pumping of the cavities is experimentally preferable. 
We especially note that since the parameter C cx N^ 1 , 
by increasing such cavity pumping C can, in principle, 
be made arbitrary small implying strong nonlinearity and 
chaotic dynamics. Consequently, by measuring the inten- 
sity of the outgoing field, the characteristics of the evolu- 
tion (Rabi, Josephson, chaotic) should be easily achiev- 
able with current state-of-the-art optomechanical experi- 
ments [llT ] . The analysis is purely semi-classical and any 
reservoir induced fluctuations have been neglected [21|. 
In the large excitation limit this cannot be a shor ting 



as verifie d by the current phonon lasing experiment 22] 
(see also 23] ). In the more purely quantum regime how- 
ever, decoherence outings may become more significant 
for the coherent tunneling process. In a full quan- 
tum analysis for a similar system was presented and it 
was found that for external pumping the Josephson ef- 
fects survives despite coupling to surrounding reservoirs. 
Moreover, most recently a photonic analog of the intrin- 
sic Josephson effect was studied and it was found that 



even deep in the quantum regime, the Josephson effect 
survives around 10 full oscillations 24J. The same is sup- 
posedly true in the present system, and we have verified 
numerically that the same type of mean-field dynamical 
structures are encountered for coherently driven cavities. 
Beyond the results of the present work, purely quantum 
effects like squeezing and anti-bunching are most cer- 
tainly interesting and will be studied in a forthcoming 
publication. 
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